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Abstract 

We report theoretical and experimental study on Auger decay of excitons in CU2O in a wide 
range of temperatures, from 2 Kelvin to 325 Kelvin, in single-photon surface excitation. We find 
that the Auger constant linearly increases with temperature above 100 Kelvin. This constant is 
important for experiments on Bose condensation of excitons. 
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I. INTRODUCTION 



There is a long history in investigating excitonic effects in naturally-grown high-quality 
CU2O, because of the inversion symmetry in CU2O and the large binding energy of the exci- 
tons, 150 meV. The inversion symmetry forbids any direct dipole recombination of excitons 
which gives excitons long lifetime up to microsecond, and 150 meV is equivalent to a tem- 
perature of 1740 Kelvin which allows excitons to exist at room temperature^. We are only 
interested in the lowest excitonic state, the "yellow" exciton series, consisting of electrons 
in the Tq level of the conduction band and holes in the Tj level of the valence band with 
2.173 eV gap energy. By electron-hole exchange, the "yellow" excitons split into a triplet 
orthoexciton and a singlet paraexciton, which lies 12 meV lower-. 

In the last ten years, there have been a lot of experiments done with excitons in CU2O, 
either with an intense surface excitation which creates a highly non-equilibrium exciton gas, 
or with an inhomogeneous strain which confines excitons in a harmonic potential trap with 
almost constant volume. One interesting topic among those experiments is the importance 
of the Auger non-radiative recombination process in Bose condensation of excitons. In 
the Auger process, two excitons collide, and one exciton ionizes, taking the energy of the 
other exciton which recombined. Some groups believe that the density- dependent Auger 
process is a barrier for Bose condensation of excitons, since it not only gives a severe limit in 
exciton density, but also heats up the exciton gas^^^^, while others claim that the process 
is negligibleS". Trauernicht et al. first gave strong evidences of an Auger process by trapping 
the excitons with an external stress^. 

Our previous work^ has shown that the Auger recombination process is important and 
increases roughly as the square of the applied stress, but it left one question in doubt, why 
there is a very small Auger effect at stresses lower than 1.5 kbar at 2 Kelvin, while a strong 
Auger effect must occur in surface excitation with a negligible stress at higher temperatures, 
for example at room temperature. Does the surface play an important role in the Auger 
effect, or is the Auger effect temperature dependent? In this work, our theoretical and 
experimental results show that the Auger decay rate increases linearly with temperature, 
which answers the above question. 
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II. EXPERIMENT 



Excitons are created at one [110] surface of a 4.9 x 3 x 3 mm naturally-grown, high-quality 
CU2O crystal with single-photon near-resonant excitation under negligible applied stress in 
the high temperature region, 115~325 Kelvin. An imaging lens collects exciton lumines- 
cence from the same [110] surface and focuses the light onto an entrance slit of a |-meter 
spectrometer which is connected with a CCD camera and a photomultiplier (PMT). The 
exciton population dynamics are examined by means of time-correlated single photon de- 
tection. This experimental approach has several advantages. First, in the high temperature 
region, paraexciton-to-orthoexciton conversion leads the excitons to reach thermal equilib- 
rium quickly. Second, as pointed out by Warren and Wolfe^., the diffusion constant at high 
temperature is much smaller than that at low temperature, which allows a constant-volume 
measurement. Also, in these experiments, we use a cavity-dumped ultra-fast dye laser which 
has a long period, 260 ns between pulses. Therefore paraexcitons definitely reach thermal 
equilibrium with orthoexcitons. The bath temperature in the experiments is controlled by 
a homemade PID controller with an accuracy of ± 2 Kelvin. 

In CU2O, there are two possible single- photon transitions, a direct transition by emitting 
a photon, and a phonon-assisted transition by emitting a photon and an optical phonon. 
Because the indirect transition is allowed for all exciton states, the indirect phonon-assisted 
luminescence spectrum shows the energy distribution of excitons^. At high temperature, 
the direct single-photon transitions for orthoexcitons and paraexcitons are almost forbidden, 
while the phonon-assisted process depends weakly on temperature-. In Section 4, we will 
talk about phonon-assisted orthoexciton recombination processes, assisted by all possible 
optical phonons and acoustic phonons, although the orthoexciton creation process via a 
Fj'g longitudinal optical phonon is much stronger than other phonon-assisted orthoexciton 
creation processes^. 

In the surface excitation, the absorption coefficient of CU2O has a strong relation to 
bath temperature and laser photon energy^ii^. We calculate absorption lengths of CU2O 
for the input laser at 606 nm and for the orthoexciton luminescence at each temperature 
from Pastrnyak's reporlj^ and J. B. Grun's paper^° on the absorption coefficient of CU2O, as 
shown in Figure 1. With a focused laser spot with radius of 120 fim, the orthoexciton has 
a volume of 7rx(120 yum)^x (Absorption-length), which is in the order of 10~^cm~^. As we 
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did in previous worksii^, we assume that one absorbed photon produces an exciton, since 
the laser is tuned close to the exciton resonance. At the highest laser power of 2.64 nJ per 
pulse, with 30% absorbed in the crystal after considering all the optics involved, there are 
10^ orthoexcitons created initially, which implies an initial orthoexciton density in the order 
of lO^^cm"^. 

Because the absorption length at the orthoexciton luminescence line is of the order of a 
millimeter at high temperatures, which is comparable to the width of CU2O sample, we have 
to add a correction factor to deduce the exciton density from the luminescence intensity at 
each temperature. The correction factor, CF, is equal to the ratio between Idetected, the real 
detected luminescence intensity, and Iq, the ideal detected luminescence intensity without 
taking into account of the absorption of the luminescence in CU2O, by a Photomultiplier. 

If we assume that one absorbed photon creates one exciton, the total number of excitons 
at a certain position x in the sample, is proportional to the number of absorbed photons. 
Tlx oc (e^i ), where < x < L and L, the width of CU2O, is 3 millimeter. Li and L2 are 
the absorption lengths for the laser and orthoexciton luminescence at each temperature, 
respectively. Because of the absorption of the exciton luminescence in the sample, the 
number of the detected exciton from the position x is proportional to (e"^). In other words, 

— X —X 

the total contribution to the luminescence from the position x is proportional to (e^i e^a ). 

Therefore, the detected luminescence intensity is proportional to the integration of the 
detected excitons over the whole length of the sample, and the correction factor is calculated 
as the following, 

CF=Lt!££^= f\^e^dx - f\^dx (1) 
Iq Jo Jo 

^CF = ^.^^. ' (2) 

^1 IT + 1 - e ^1 

In principle, we should account for re-emission from the reabsorbed photons, but since 
the re-emission mostly goes in directions other than x (see inset of Figure 1), this will be a 
small correction. 

III. RESULTS AND DISCUSSIONS 

Figure 2 shows the typical orthoexciton photoluminescence intensity as a function of time 
after a short laser pulse, at two different temperatures, 150 and 325 Kelvin, and under three 
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different pulse energies, 2.64 nJ, 2.64/2=1.32 nJ and 2.64/6=0.44 nJ, multiplied by overall 
factors of 1, 2 and 6 respectively. At each temperature, the different decay times for different 
laser powers indicate that there is the density-dependent Auger process of excitons. 

Our previous worki*^ presented coupled rate equations including the Auger decay rate, 
which explained well the data at 2 Kelvin under different stresses and laser pulse energies. 
However, because of phonon emission and absorption, the orthoexciton and paraexciton are 
well coupled in the high temperature region^, and therefore the coupled rate equations can 
be modified to the following: 



In these equations, indies o and p stand for "orthoexciton" and "paraexciton", respectively, 
A is the Auger decay constant, n is the exciton density, r is the decay time or conversion 
time, and the indies o-p and p-o mean orthoexciton-paraexciton conversion and paraexciton- 
orthoexciton conversion, respectively. In each Auger process, two excitons collide and end 
up with one exciton recombining and the other exciton ionizing. Therefore, only the ionized 
exciton, half of the excitons participating in the Auger process, will be added back to exciton 
population. Because the orthoexciton is a triplet state and the paraexciton is a singlet state 
and spin is randomly selected in ionization, | of the ionized excitons in the Auger process 
will be returned as orthoexcitons and | of these excitons will be returned as paraexcitons. 
We have to consider inter-conversions between paraexcitons and orthoexcitons in the high 
temperature region, although only orthoexcitons are converted to paraexcitons at 2 Kelvin^. 
Because energy splitting between orthoexciton and paraexciton at the zone center is 12 meV, 
and orthoexciton and paraexciton reach thermal equilibrium rapidly at high temperature^, 
the conversion time in the inter-conversion mechanisms has a relationship. 



All the data, under ten different temperatures and three different laser pulse energies at 
each temperature, are well fit with the above coupled equations, as shown on the solid lines 
of Figure 2 (a) and (b). The fit results imply that the Auger constant of the orthoexciton is 
equal to that of the paraexciton at each temperature. The results also show that the radiative 
lifetimes of paraexcitons and orthoexcitons are on the order of several hundred nanoseconds, 





(3) 



= r, 




(5) 
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and the inter-conversion time between paraexcitons and orthoexcitons is on the order of 0.1 
nanosecond which is consistent with the report from Wolfe's group^ and our previous results^. 
Because paraexcitons rapidly reach thermal equilibrium with orthoexcitons, in the fitting 
process we can assume that at all times after t = 0.1 ns, 

12meV 

np{t) = Uoit) xe'^BT . (6) 

An interesting result from the fitting is that the Auger constant of excitons linearly 
increases with temperature in the region from 115 Kelvin to 325 Kelvin, as shown in Figure 
3. Several groups have studied the temperature dependence of the Auger decay process 
in CU2O, but the highest temperature they used is 77 Kelviii^iii. Their works show that 
the Auger constant of orthoexcitons is almost independent of temperature. The above two 
results are not inconsistent, which is explained in detail in next section. 



IV. PHONON-ASSISTED MECHANISM 

Kavoulakis and BaynJ^ introduced a model based on Fermi's Golden rule and k-p theory 
to understand the Auger recombination process with an optical phonon-assisted mechanism 
at low temperature. We follow their steps, but modify it with considerations of high tem- 
perature effects and all possible phonons in an Auger process, in which two excitons with 
momentum K and P collide, producing one ionized electron-hole pair with momenta k^. and 
khi and absorbing or emitting a phonon. 

In the phonon-assisted process, ^y^"' , the Auger decay rate per unit volume is a func- 
tion of Tk,p and fx, fp, the decay rate and distribution functions of the two excitons, 
respectively: 

^F = V^fKfp^^,P-An^ , (7) 

^ ^ K,P 

Where A is Auger constant for the exction and n is the density of the excitons. 

From Fermi's Golden Rule, the decay rate of two excitons with momenta K and P is 

2n 

X [{uph + l)5iEK + Ep- hujQ - ^cfc, - i^v,kj 

+nphS{EK + Ep + hujQ - ec,k, - e^^kj] , (8) 
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The first term in Equation 8 is for a process emitting a plionon witli momentum Q, wliile 
tlie second term is for a phonon absorption process by absorbing a phonon with momentum 
Q. The term = Ij I^c^p^/'^^bT) _ jg ^j^^ density of the phonon at temperature T, and 
fcs is Bohzmann's constant. 

In the equations, based on k-p theory^^, the matrix element of an optical phonon-assisted 
process, yiopUcah is given by the following equation, 

^Vloptical ~ ^ TT X X — — 

x( "^''"^ )^ -x5^ (9) 

where T>opticai is the deformation potential for an optical phonon, EopUcai is the energy of 
the optical phonon, Ec and Sc' are energies of the lowest conduction band and the only 
odd parity conduction band Fg , respectively, is the energy of the valence band Fy , P^^c' 
is the matrix element of the momentum operator between the conduction and valence band 
at zone center, p is the mass density of CU2O, eoo is the dielectric constant, and a^ is the 
exciton Bohr radius. 

In an acoustic phonon assisted process, kmcessary the necessary momentum for an acoustic 
phonon is to satisfy with energy and momentum conservations between the lowest conduction 
band and the only odd parity conduction band, that is to say, kmcessary = (^c — £c')/{ft-^) 
= 1.515x10^ cm^^. In the discussion, u is the acoustic phonon speed. However, kmax, the 
maximum momentum of an acoustic phonon, is equal to ^ = 7.375x10^ cm~^ because 
it is restricted to be in the first Brillouin zone, where is the lattice constant of CU2O. 
Therefore there is no possible process assisted by the acoustic phonon and from now on all 
discussions will only consider the process assisted by the optical phonon. 

As mentioned by Kavoulakis^^, Ek, Ek^, and Ek^, the energy of an exciton with momen- 
tum K, the energy of an electron in the conduction band, and the energy of a hole in the 
valence band, respectively, are given hy^, 

Ek = E,,p + -Aj^ (10) 

^' '''exciton 

h^k^ 

E,^ = E,^p + ^ (11) 
Ek, = ^ (12) 
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Because the momenta of excitons, K and P, are determined by exciton thermal motion, we 
can make a reasonable assumptionii, 

K,P<^K,kh . (13) 

By substituting Equation (9) through (13) into Equation (8), Tk,p the decay rate of two 
excitons in the phonon-assisted process is equal to 

9I4 3p4 t r)2 IP 12 

Yphonon-assisted z, /I e /X/i ^ ^optical ^ \-'^c',v\ 

^'^ el^a\V IV? E optical P {^0,0 - ^dflY {^dfi - ev,oY 

[1 + {Kaeff e^o,u.ai/kBT _i\ ' ^^^) 

The integration is restricted to the first Brillouin zone. All constant parameters used 
in the calculation are listed in Table 1. The optical deformation potentials and |Pc,i)P are 
estimated from fitting of experimental data. 

Therefore, by substituting Equation (14) into Equation (7), we get Aphonon-assisted, the 
phonon-assisted Auger constant, in the unit of cm^/ns, 

Aphonon-assisted ^ 3.698 X 10"^^ X [l + ^^^J^ _ J + 1.298 X 10"^^ X [l + J . (15) 

Only two optical phonons, rj"2 and T^^, make important contributions. 

We plot Equation 15, a theoretical prediction of the exciton Auger constant assisted 
by all possible phonons from 2 to 325 Kelvin, and our experimental Auger recombination 
constants in the high temperature region ( 115 ~ 325 Kelvin ) in Figure 3, which shows 
that our experimental data in the high temperature region can be well fit by the theoretical 
prediction of the Auger constant from the phonon-assisted mechanism. 

At low temperatures, where the phonon energy is much smaller than the thermal energy 
k^T, Uph the density of the phonon is much less than 1. Therefore only the phonon emission 
process is important and phonon absorption process is negligible. We can approximate: 

^K,P = X E |M2j(l-n,,,J(l-n„,,J 

X 6{Ek + Ep - hujQ - ec,k, - ev,ki,) (Low T) (16) 

After substituting Equation (16) into Equation (7), we find that the Auger decay constant 
is weakly dependent on temperature when the temperature is lower than 77 Kelvin, which 
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is shown in Fig. 3 and also consistent with the reports from Warren et al^ and Kavoulakis 
et alM^. 

Therefore, we can exphcitly state our conclusion on the temperature dependence of the 
Auger recombination constant of excitons as the following: the Auger constant is weakly 
dependent on temperature when the temperature is below 77 Kelvin, and linearly increases 
with temperature in the high temperature region. 

To verify our model, we make two comparisons. The first comparison is between a 
reported Auger constant from Wolfe's group^ and our theoretical prediction at 77 Kelvin, 
which are consistent with each other, though our theoretical prediction is about 4 percent 
larger, as shown in Figure 4. The other comparison is made between the experimental 
data and our theoretical model at 2 Kelvin. Our theoretical result indicates that the Auger 
constant does not continuously decrease to zero at zero stress at 2 Kelvin, but becomes 
a constant, 4.996 x 10~^^ cm^/ns, which is comparable to the experimental data at 1.5 
kbar at the same temperature reported in the previous work^, as shown in Fig. 4 also. 
The comparison implies several things. First, the Auger constant has a negligible stress 
dependence when the stress is lower than 1.5 kbar. Second, the two results are consistent, 
although the theoretical prediction from our phonon-assisted mechanism is about 9 percent 
larger than the extrapolated data. Third, the surface doesn't play an important role in 
the Auger effect. Since the 4 to 9 percent discrepancy between experimental data and 
theoretical prediction is in the acceptable experimental error range, our modified phonon- 
assisted mechanism may be a good model for the Auger recombination rate of excitons in 
both low and high temperature regions. 

V. CONCLUSIONS 

From the above theoretical and experimental studies, we may say that in surface 
excitations, the Auger recombination constant linearly increases with temperature at high 
temperatures, and weakly depends on temperature in the low temperature region, which is 
consistent with our modified model of phonon-assisted Auger recombination. However, we 
do not have a model for the stress dependence. 
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Parameter 


Value 


References 


^C30 


6.46 


11 


p 


6 gram/cm^ 


13 


V 


4.5x10^ cm/second 


12 


me 


ml 


4 


rrih 


0.7m° 


4 


rrig 


9.1x10^* gram 


14 


as 


7x10""^ cm 


10 


aL 


4.26x10-^ cm 


10 




l.lOiO.Ol eV/A 


this work 


i 15 


0.75±0.01 cV/A 


this work 


\Pc',v\/n 


0.55±0.01 A"^ 


this work 




0.449 eV 


11 


Ed 


2.622 eV 


11 




13.8 meV 


15 




18.7 meV 


15 



TABLE I: Constant parameters used in the Auger decay rate calculation. 
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FIG. 1: Absorption lengths of the input laser at 606 nm ( black squares ), and the exciton 

luminescence ( open triangles ) in CU2O at each temperature in the high temperature region. 
Inset: Schematic of the experimental set-up in this work. 
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FIG. 2: (a) Orthoexciton luminescence decay at T = 150 Kelvin for three different pulse energies, 
2.64 nJ, 1.32 nJ and 0.44 nJ, as shown in the figure. Dots: Phonon-assisted orthoexciton lumines- 
cence intensity as a function of time after a short laser pulse. Solid lines: fit to the coupled rate 
equations discussed in the text, (b) Orthoexciton luminescence decay, as in (a), but at T = 325 
Kelvin. 
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FIG. 3: Predicted and experimental Auger constants of the excitons as a function of temperature 
in the high and low temperature regions. 
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FIG. 4: Stress and temperature dependence of the orthoexciton Auger constant. Open triangles: 
data reported by Snoke et al} at 2 Kelvin for three different stresses, 1.5, 2.5 and 3.5 kbar. Open 
stars: data reported by Wolfe et al? with negligible stress at 77 Kelvin. Black triangles and black 
stars are predictions from our theoretical model with negligible stress at 2 Kelvin and 77 Kelvin, 
respectively. 
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